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Abstract. We consider a model in which massive stars form in a self-gravitating accretion disk 
around an active galactic nucleus. These stars may evolve and collapse to form compact objects 
on a time scale shorter than the accretion time, thus producing an important family of sources for 
LISA. Assuming the compact object formation/inspiral rate is proportional to the steady-state gas 
accretion rate, we use the observed extra-galactic X-ray luminosity function to estimate expected 
event rates and signal strengths. We find that these sources will produce a continuous low-frequency 
background detectable by LISA if more than ;> 1 % of the accreted matter is in the form of compact 
objects. For compact objects with m ^ IOM0, the last stages of the inspiral events should be 
resolvable above a few mHz, at a rate of ~ 10—100 per year 

INTRODUCTION 

One of the major challenges that will face gravitational wave (GW) astronomers is the 
successful discrimination between instrumental noise, stochastic GW backgrounds, and 
individual, resolvable GW sources. For example, the population of galactic white dwarf 
binaries in close orbits will provide a major contribution to the LISA noise curve in 
the 0.1 — 1 mHz band JU]. At the same time, this confusion "noise" can also be treated 
as a signal, and its shape and amplitude will provide important information about the 
distribution and properties of white dwarf binaries in the galaxy. 

It has recently been proposed that, in the self-gravitating accretion disks of active 
galactic nuclei (AGN), massive stars could form and evolve, eventually collapsing into 
compact objects and merging with the central black hole [2, 3]. At different stages in 
the inspiral evolution, this population will contribute to the GW background confusion 
or alternatively produce individual, resolvable chirp signals. It is therefore a matter of 
theoretical and practical interest to understand the nature of such a population. 

In this paper we derive a relationship between the observable electro-magnetic (EM) 
emission and the predicted GW emission from AGN, following the procedures and 
summarizing the results outlined in our longer paper 0]. In particular, we use the hard 
X-ray luminosity function of Ueda et al. |0] to infer the accretion history of supermassive 
black holes (SMBHs) out to redshifts of z ~ 3. Then we assume a few simple scaling 
factors, such as the average (Eddington-scaled) accretion rate and the efficiency of 
converting accretion energy to X-rays, and derive the GW spectrum that might be seen 
by LISA. 

Depending on the specific model parameters, we find this background could be an 



TABLE 1. Glossary of dimensionless parameters, with allowable and preferred values 



symbol 


min 


max 


preferred 


description 


/acc 





1 


1 


fraction of SMBH mass due to accreted gas 


,/co 





1 


0.01 


fraction of SMBH mass due to accreted compact objects 


fx 





1 


0.03 


fraction of EM radiation in X-rays 


/Edd 





^ 1 


0.1 


typical fraction of Eddington luminosity/accretion rate 







1 


0.2 


efficiency of converting accreting gas to EM radiation 







1 


0.2 


efficiency of converting compact objects to GW radiation 



important class of LISA sources, similar in strength and event rates to extreme mass-ratio 
inspirals from captured compact objects [6]. As in those sources, here too it is a matter 
of preference as to whether the background should be thought of as signal or noise. But 
for higher masses (perhaps as large as m ~ IO^Mq [|2lfl), disk-embedded compact objects 
should produce individual, resolvable inspiral events with high signal-to-noise over a 
wide band of frequencies. 



THE HARD X-RAY LUMINOSITY FUNCTION OF AGN 

We begin with a short discussion of notation. The results derived below include a number 
of dimensionless parameters, most of which can take values between and 1. We 
divide these parameters into two general classes: efficiencies and fractions. Efficiencies, 
denoted by 7], are believed to be determined by more fundamental physics, and typically 
have more stringent lower- and upper- limits. Fractions, denoted by /, are more model- 
dependent parameters and less-well known than the efficiency parameters, and thus have 
a larger range of acceptable values. A summary of these model parameters appears in 
TableUl 

A growing consensus has been forming that SMBHs grow almost exclusively by 
accretion, suggesting /acc — 1 (see, e.g. A corollary of this assumption is that most 
AGN should be rapidly spinning, with dimensionless spin parameters of a/M ~ 0.9 — 
0.998, giving EM efficiencies for a radiative disk of r^em — 0.15 — 0.32 [8]. Similarly, 
assuming the compact objects are on circular, adiabatic orbits, we set r/gw = Tjem, 
ignoring any EM or GW emission from the plunging region. 

We write the X-ray luminosity Lx as a fraction fx of the bolometric luminosity, which 
in turn is a fraction /Edd of the Eddington luminosity L^dd- 

Lx = fxfEddLEdd{M) = fxVemM^ccc'^. (1) 

The Eddington limit is a function only of the SMBH mass: L^ddi^) = 1.3 x 
10^^{M/Mq). Over the range of redshifts and luminosities we are probing, typical 
accretion rates are estimated to be /Edd ~ 0.1, but could conceivably be even greater 
than unity L2J- Since the energy density in the cosmic infrared background (also 
dominated by AGN) is about 30 times greater than the X-ray background, we set 
fx = 0.03 \ lu\. Lastly, the parameter /co is the fraction of total mass accreted in the 
form of compact objects. Since the astrophysical mechanisms that actually determine 



this fraction are not yet well understood, we set it to a conservative value of 0.01. If it 
were much higher, the disk would be entirely fragmented and thus not efficiently emit 
EM radiation. And as we will see below, a value much below 0.01 would result in a GW 
signal undetectable by LISA. 

We approximate the intrinsic (i.e. directly produced by the accretion process, and 
before reprocessing and/or partial absorption within the host galaxy) X-ray luminosity 
function per comoving volume after Ref. Js]], roughly giving a broken power-law dis- 
tribution with higher average luminosities at higher redshift. The fiducial values for the 
model parameters and the observed range of luminosities correspond to SMBH masses 
in the range 1O^M0 ^ M ^ IO^'^Mq, consistent with the masses inferred from observa- 
tions of velocity dispersions. 



THE GRAVITATIONAL WAVE SPECTRUM 

We begin by considering the inspiral of a single compact object of mass m onto a SMBH 
of mass M ^ m and specific angular momentum a. Using geometrized units such that 
G = c = 1, a particle on a circular, equatorial orbit around a Kerr black hole has an 
orbital frequency (as measured by an observer at infinity) of 



and specific energy 



m 



r{r^ - 3Mr + lay/Mr) 1/2 ' 



Thus the total energy emitted in gravitational waves down to a radius r is Eg^{r) = 
m — E{r). The GW energy emitted between frequency / and f + df for such an event is 



dEg^ji^ dEgyi^ f d f 



df dr \ dr 



(4) 



Here, we restrict the GW emission to the quadrupole formula for circular geodesic orbits, 
thus we consider only GW frequencies at twice the orbital frequencies (/ = 2/orb)- 

We will generally want to restrict equation Q to a range of frequencies /min < / < 
/max, where /mm is determined by the LISA sensitivity and /max is the GW frequency at 
the inner-most stable circular orbit (ISCO). The ISCO frequency in turn is determined 
solely by the SMBH mass and spin, giving /max ~ 4 — 30 mHz for M = \O^Mq. In 
Figure [l] we show the characteristic strain spectrum for a single inspiral event for a 
range of black hole masses and spins. 

Integrated over redshift z, the observed GW energy density per logarithmic frequency 
is given by (e.g. [11.1 ') 
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FIGURE 1. Characteristic GW strain amplitudes for individual inspiral events, where a black hole with 
m=l OMq merges with a SMBH of mass M = 1 0^ , 1 0^ , 1 O^Mq at a redshift of z = 1 . For each value of M, 
we show the spectra for two spin values, a/M — (solid) and 0.95 (dashed). The dot-dashed and dotted 
lines are the sky-averaged LISA noise curves with and without the contributions from galactic binaries, 
respectively. 



where /, = (1 +z)f, and R{z) is the rate of inspiral events per comoving volume. 
Cosmology enters with the term \dt /dz\ = [(1 where for a flat geometry, 

H{z) = Ho [auii+zf + a^] ^'^ . (6) 

Throughout this paper we will assume a standard ACDM universe with = 0-3, 
= 0.7, and Hq = 72 km s^^ Mpc^^. 
The event rate for a single AGN is simply the X-ray luminosity divided by the total 
X-ray energy emitted between inspiral events: Ex = £^aw/x//co- Integrating over the 
luminosity distribution function, we get 

fx J dlnLx 

where dn{Lx,z) /dlnLx is the intrinsic luminosity distribution function in units of 
Mpc~^. Combining equations ^ and Q, the total (time-averaged) gravitational wave 
spectrum is 

dLx^—^ia (8) 
JlnLx £'gw dlnfz 

The GW spectrum Eg^{f) from each individual AGN is a function of the SMBH mass, 
which in turn is determined by the X-ray luminosity through equation ([T]). Note that the 
integrated spectrum is independent of m, as long as m is small enough so that the inspiral 
waveform cannot be individually resolved. One measure of this resolvability is the duty 
cycle, described in the next section. 
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FIGURE 2. Dependence of the GW background on the model parameters /edd (left) and /co (right). 
The fiducial values for the model give a background signal comparable to the LISA noise curve, including 
unresolved galactic white dwarf binaries. 



Following Refs. 0,112], we will want to compare directly the background defined 
in equation (HJ) to the spectral density of the detector noise Sn{f), which has units of 
inverse frequency. In this case, ^/fS^Jj) will be a dimensionless strain. Averaging over 
the entire sky, weighted by the LISA antenna pattern, gives 



4 1 Jpgw(/) 
np dlnf 



W) = ^^^^^- (9) 



Throughout the paper we use the so-called sky and detector averaged instrumental 
spectral density for LISA, augmented by the white-dwarf galactic confusion noise, as 
given in [6]. 

In Figure |2l we plot the nominal GW background from equation (HJ) with the model 
parameters listed in Table 1. Also shown are the effects of varying the accretion rate 
/Edd (Fig. Eh) and the fraction of accreted mass in compact objects /co (Fig. |2b)- By 
increasing /Edd, the effect is to reduce the AGN mass inferred from equation dU) and 
thus increase the GW frequency of the signal, shifting the curves to the right. From the 
leading term in equation (HI), it is clear that the total GW power is simply proportional 
to /co, so increasing this parameter linearly increases the total amplitude of the GW 
spectrum. 



EVENT RATES AND RESOLVABLE SIGNALS 

The GW spectra calculated from equation (EJ) represent the time-averaged signal from 
all the AGN in the universe, but at any single time, there may only be a few sources that 
are emitting at a given frequency. This can be seen from comparing Figures [l] and |21 
while the signal-to-noise ratio of a single inspiral can be very large around 1 — 10 mHz, 




a relatively small fraction of the inspiral time is spent in that band, significantly reducing 
the time-averaged strain amplitude. 

One way of estimating the number of individual signals at a given frequency is by 
calculating the duty cycle £>(/): 

D{f) = / ^^T^ / dLx-f-tUa (10) 

where the cosmological volume element is dV /dz = ATtr^{z) /H{z) and the "coherence 
time" ?nnh( f) is approximated by the Newtonian limit of the radiation reaction formula 
(e.g. iH): 

In Figure El we show the duty cycle for the fiducial model parameters and a range of 
compact object masses m. The sharp cutoff around 40 mHz is not physical, but rather 
due to the somewhat artificial low-end cutoff in the luminosity function, corresponding 
to a minimum value for M and thus maximum attainable frequency. Note that D{f) 
is proportional to m^^, since smaller m means more compact objects, and also slower 
inspiral rates, thus each source spends more time around a given /. 

It is important to understand that D{f) represents the total number of sources in the 
observable universe, not all of which would be individually resolvable with LISA. To 
estimate the subtractable portion of the signal, we integrate the total signal-to-noise ratio 
(SNR) for a three-year mission lifetime (see, e.g., [12]). For a given m, M, and source 
distance d, we can calculate the minimum frequency above which the source would be 
detectable with SNR above a certain threshold (here we use 15). Any contribution to 
equation dSj) above that frequency is ignored, resulting in an unresolvable GW spectrum. 
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FIGURE 4. The fiducial time-averaged signal (top curve), along with the unresolvable portions of the 
GW spectra, after subtracting out individual events with integrated SNR above a threshold of 15. 

plotted in Figure IH As expected, with increasing m, more inspiral events are detectable 
and the remaining background confusion noise is diminished. 
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